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1. Introduction 



One Matrix Models are quantum mechanical systems whose partition function is de- 
fined by an integral of the form: 

f N 

Z M = J dM exp{--trF(M)} , (1) 

where M is an iV x N matrix and the potential V(M) is a polynomial in M. In the last few 
years, there has been tremendous progress []l]-fi in matrix models through the discovery of 
a connection of a certain class of these models to two dimensional gravity coupled to (p, q) 
minimal confer mal matter. This happens when M is a hermitian matrix (HMM) or when 
one considers generalizations to a (q — 1) hermitian multi-matrix model (MHMM), where 
(q — 1) hermitian matrices are coupled linearly to each other 0. In the former simpler 
case, the Feynman graphs of the zero dimensional field theory are viewed as being dual to 
a discrete dynamical polygonation of an oriented two dimensional Riemann surface. Then 
the perturbation series can be summed in the form 



oo 

= Y.mz h (2) 



h=0 

where X = ^ — E + L is the Euler character of the corresponding surface and h is its genus 
given by x = 2 — 2h. Since the number of vertices V, edges E and loops L of the Feynman 
graph correspond respectively to the number of faces F, edges E and vertices V of the 
dual graph, the above series can be shown to correspond to the discretized version of the 
partition function of pure two dimensional gravity 

In (|D A is the area of the surface, [ib and Gb the bare cosmological and Newton's constant 

and C(T) is the symmetry factor of the polygonation corresponding to dividing by the 

volume of the diffeomorphism group of the surface. The equality of @ and (|3|) is achieved 

i 

by identifying A = e _AtB and iV = e 47rG s . The action in (|3p can also be viewed as the 
action of a string theory embedded in zero dimensional spacetime 



S str = log K B J d 2 ^R + /j B J d 2 i^. (4) 



Then (|2]) gives the genus perturbation expansion with Kb = jr, the bare string coupling. 
The naive continuum theory is taken by letting iV — > oo. In this case the area A diverges 
and the polygonated surface is thought to approach a smooth Riemann surface. For a 
critical value \x c of [xb the increasing entropy of large surfaces compensates the Boltzmann 
factor and the system undergoes a (third order) phase transition. If the critical point is 
approached in an arbitrary way, only the sphere Zq contributes to (0). The remarkable 
observation [l]-|4|] was that since the singular part of Zh ~ (^b — /U c ) x< - 1+ '^* ) with k a 
positive integer, one can obtain contributions from all genera by simultaneously taking 
the large iV-limit and letting \xb approach its critical value \i c in a coordinated way. The 
integer k labels a series of multicritical points reached by tuning k parameters in the 
potential V(M). Introducing a cutoff a in the theory, we define the string coupling kq and 
renormalized cosmological constant (ir to be 

«o = N , m = q2 ■ (5) 

The double scaling limit is defined by taking N — > oo and [ib — ► At c while keeping k and 
\ir fixed. Then the continuum limit of (|!|) becomes 



j> x Z fc , (6) 



^str — 

h=0 

with k = The series (^) is horribly divergent. It is non-Borel summable since every 
term increases as (2h)l. This reflects our ignorance in summing the perturbation series of 
string theory although the fixed genus partition function Zh can be calculated and is well 
defined. Happily, the theory is exactly solvable at the multicritical points and its dynamical 
content is given by a single differential equation, the string equation. The string equation 
is a differential equation in the variable x satisfied by the specific heat — d 2 log Z, with 
h? = x~( 2+ ^\ It possesses solutions that in the weak coupling limit k — > are asymptotic 
to (|D and we say that the double scaling limit provides a non-perturbative definition of 
Z s tr- Indeed comparison with calculations directly from the continuum theory indicates 
that Z s tr corresponds to two dimensional gravity coupled to (2k — 1,2) minimal conformal 
matter. Even more interesting is the discovery that the double scaling limit of (q — 1) 
MHMM gives two dimensional gravity coupled to (p, q) minimal conformal matter || . 

Unitary One Matrix Models (UMM) form another interesting class of matrix models. 
These are defined by (|1|) with M being a unitary matrix U. The interest in those models 



arose a long time ago when Gross and Witten [BJ showed that the partition function of 

V 

two dimensional U{N) QCD on a lattice is given by Zqcd = {Zu)~°? an d that the theory 
undergoes a third order phase transition in the large iV limit (V is the volume of the two 
dimensional world and a is the lattice cutoff). The theory was also shown to posses a 

2k \ 

A AT Ol.J.1 (-.n^ n i / I " 



2k + l 



double scaling limit N — > oo and A — > A c with t = (1 - -^iV 5 ^ and y = (1 - x")iV 
held fixed |7]J§. The string equation is a 2k th order differential equation of the function 
v in the variable x = t + y, with v 2 = —d 2 log Z. It has solutions that are asymptotic to 
© in the limit x — > oo with k 2 = £ _< - 2+ r). The identifications of those solutions with 
conformal field theories coupled to two dimensional gravity or other interesting systems 
is still, however, an interesting open problem. Some interesting suggestions have been 
made in ||. Moreover, the surface interpretation of UMM is not as clear as in the case of 



HMM. In [[10] Neuberger views the unitary matrix as U = e lM where M is hermitian and 



introduces N x N hermitian fermionic matrices ip and ip to exponentiate the Haar measure 
dU — ► dMdet(jjf). The resulting surfaces contain an infinite number of types of bosonic 
vertices forming bosonic "webs" and fermionic loops forming their boundaries that might 
allow a stringy interpretation of the UMM. For another interesting suggestion see 1TTJ] . It 



is also interesting to note that UMM belong to the same universality class as the HMM 
in a different class of multicritical points, the double-cut HMM [PJI^j. This is expected 
since the critical behaviour is governed by the scaling of the density of the eigenvalues at 



the edge of its support jl3| and the eigenvalues of the two models models scale identically 
there. 

The continuum theory obtained in the double scaling limit has a very rich mathemat- 
ical structure. When one considers perturbations by the scaling operators < > with 
sources t^, the dependence of the specific heat (or its square root for the UMM) on the 
"times" tk is given by KdV flows |5|,|nj for the HMM and mKdV flows for the (symmetric) 
UMM. The partition functions of the theory are found to be given by the corresponding 



r-functions |15|-[17j which can be thought as sections of a line bundle over the Universal 



Grassmannian. Furthermore the r-functions that solve the string equation are annihilated 



by constraints which for the one matrix model are the Virasoro constraints fi5|- |Iql . All 
of those results have counterparts in the discrete theory. The integrable flows are now 
with respect to the couplings in the potential V(M). For the UMM these are given by 



Toda flows on the half line |T9| and the partition function is given by the product of two 



Toda-chain r-functions. The Virasoro constraints L n have the simple interpretation of 



corresponding to invariance of the partition function under specific transformations, which 
for the UMM are given by SU = e n (U n+1 - U 1 '" 1 ). 

An interesting observation is that the string equation can be written in the form 
[P, Q] = 1 where P and Q are differential operators for the HMM || and 2x2 matrices of 
differential operators of specific order for the UMM |2(J . They correspond to the continuum 
limits of operators acting on the space of orthonormal functions used to solve the model. 
One can use this form of the string equation to determine easily the points in the Universal 
Grassmannian that solve the string equation ||21|1 . For the UMM these are found to 
correspond to a pair of points V\ and V2 in the (big cell of the) Sato Grassmannian 
satisfying certain invariance conditions. It is very important that the mKdV evolution of 
V\ and V2 gives new solutions to the string equation. The r-functions that correspond to 
V\ and V2 are shown to satisfy the Virasoro constraints in this formalism since the 
constraints are derived from the same invariance conditions that solutions to the string 
equation satisfy [^-^(3 . 



This article is organized as follows. In section 2 we review the discrete formulation 
of the symmetric UMM. The method of the orthogonal polynomials in the trigonometric 
basis is summarized and the Toda flows and Virasoro constraints are discussed. In section 
3 we describe the double scaling limit and describe how the mKdV flows arise. In section 4 
we give a non-rigorous approach to the connection between the Sato Grassmannian and the 
mKdV flows starting from the finite dimensional Grassmannians. In section 5 we describe 
the connection of the Sato Grassmannian to the solutions to the string equation. 

2. The Symmetric Unitary Matrix Model 

In this paper we will study the UMM defined by the one matrix integral 

/TV 
DUexp{--TrV(U + U^)}, (7) 

where U is a 2N x 2N or a (2N + 1) x (2N + 1) unitary matrix, DU is the Haar measure 
for the unitary group and the potential 

V(U) = J2dkU\ (8) 

A;>0 

is a polynomial in U. As standard we first reduce the above integral to an integral over 
the eigenvalues |||27| Zj of U which lie on the unit circle in the complex z plane. 



z *= /{n^}i A ^i 2ex p{-xE^+<)}' ( 9 ) 
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where A(z) = Y[ ( z k ~ z j) is the Vandermonde determinant. The Vandermonde determi- 

k<j 



nant is conveniently expressed in terms of trigonometric orthogonal polynomials [28 



z n ±z~ 



a n,n-A Z 



± Z 



-n-\-i s 



1=1 



(10) 



= ±ct(z~ 1 ) 



where for U(2N + 1) n is a non-negative integer and i m ax = n and for U(2N) n is a positive 
half-integer and i m ax = n — \ . The polynomials c~(z) are orthogonal with respect to the 
inner product 



(c + c + ) 



/ exp{-^ V(z +z*)} ct{zYcUz) 
J Lniz A 



= 5 



'ID 



(Cn,C~) =0. 

The expression for the Vandermonde determinant is 



|A(*) 



det 



where j = 1, . . . , 2iV, % 



I 3 

2' 2' 



iV - i for U(2N) and j = 1, . . . , 2N + 1, i = 0, 1, 



(12) 



AT 



for £7(2iV-|- 1) (where the line c (z) = is understood to be omitted). Then the partition 
function of the model is given by the product of the norms of the orthogonal polynomials 



J N 



n 



r (+) r (-) 

N ' N 



(13) 



The functions Tn and rA ' are Toda chain r-functions on the half line IT3 



.(-) 







2a± 



dgl 



(14) 



r (±) 

Tl + l 

JIT 



with solutions e^™ 

The orthogonal basis of polynomials chosen is especially useful for constructing the 
operator formalism of the theory. When acting on the basis of orthonormal functions 



TT^(z) 



c„(z 



(15) 
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such that 

(ir+(z).ir±(z)) = I 

2niz 

the operators z± = z ± - and z<9 z give finite term recursion relations 



(16) 



'07 



zd z T7±(z) = > ](vf) n , n+r ir^ +r (z) + {n\l^- ^r(vt) n , n }ir*(z 



2 

AT 



5^(^)n,n+r^ +r (z) + I U 



fc 



2A 



where i£ = e^"^-, Q± = e^- r± = g± = ^lZ^±ffilZ^ > , and 



Then the discrete string equation is given by the relation [zd z , z±] = z T . 
Invariance of the partition function under the transformations 

8 n U = e(U n+1 - U l ~ n ) n > 1 , 

implies that the partition function is annihilated by the Virasoro constraints 



°° d 1 d 2 

^n = Y J kg k - + - 7T^ • ( 18 ) 

^ dg k +n 2 x ^ <n og k og n -k 



k=0 



In |T^] it was argued that the string equation can be viewed as a consistency condition of 



the integrable hierarchy and the Virasoro constraints. 

3. The Double Scaling Limit 
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The continuum limit of (0) is taken by letting N — > oo. Then the eigenvalues a^, 
where Zi = e lQi , become continuously distributed over the unit circle \z\ = 1 and their 
distribution is described by the density of eigenvalues 

aa (19) 



p(a)da = 1 < a c < n . 

c 

If p(a) is given, quantities of physical interest, like the free energy F sp h = log Z, can 
be calculated. For example the saddle point approximation of (0) gives 

2 f ac f ac a — 3 

Fsph = ~r / dap(a)V(2 cos a) + P / dad3p(a)p(3) log | sin — - — | + const. , (20) 

* J-a c J-a c 2 

where P denotes the principal value of the integral. Then one can think of the eigenvalues 
as a Dyson gas of electric charges on the unit circle subject to their mutual Coulomb 
repulsion and an external potential V . In the weak coupling limit A — ► oo the eigenvalues 
tend to distribute uniformly on the circle, whereas in the strong coupling limit A — > the 
charges are localized, say at the point z = 1. The system undergoes a phase transition 
precisely when the eigenvalue distribution develops a cut at z = — 1 and it happens when 
A c = 1. Near the cut p(a) scales as 

p k (a) ~ c fc (l -sin 2 |) fc ol > ty , (21) 

and we obtain a third order phase transition with F ~ (A — A c ) 2+ ^ The k th multicritical 
point is obtained by tuning k couplings in the potential V{U) to their critical values. 

The double scaling limit [[7],|8| corresponding to the k th multicritical point is defined 
by iV ^ oo and A — > A c , with t = (1 - §)NW : y = (1 - ±)Nm held fixed. It was 
shown in [BD[ that the operators z± and zd z have a smooth continuum limit given by 

z + ->• 2 + N~W Q + z-->-2N~wQ_^ 

(22) 

zd z -> A^wt V k , 
where Q± are given by 

Q- = 

Q+- 

= a 



d + v 
d-v 

(d + v)(d~v) ^ (23) 

(d-v)(d + v) 
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and V k by 

(24) 

Here d = ^ and x = t + y. The scaling function f 2 is proportional to the specific heat 
— <9 2 InZ of the model. The operators are differential operators of order 2k. The same 
assertions hold if we introduce sources t<2k+i(ti = x) and deform the k th multicritical 

2(fc-i) 

potential Vfe to Vk(z) — J2^2l+iVi(z)N 2fc + 1 . From [zd zl Z-] = z + it follows that 

l 

[V, Q_] = 1 , (25) 

where P = — {21 + l)t 2 ;+iPz — a: with P; = P; + x. The function v(x) becomes a 
i>i 

function of x and the times {^2Z+i } and obeys the string equation 

^(2/ + l)t 2Z+1 £>R z [u] = -vx . (26) 
i>i 

where V = d + 2v, u = v 2 —v', and Rk[u] are the Gel'fand-Dikii potentials defined through 
the recursion relation 

dR k+1 [u}= (^d 3 -^(du + ud)^R k [u}, Ro[u} = ^. (27) 

The dependence of v on the times {^2Z+i } is given by the mKdV flows 

Cr -df>R k [u}. (28) 



dt2 k +i 

It is very important that (j28D is compatible with the string equation. It can be shown (see 
also section 5) that solutions to the string equation flow with (^8j) to other solutions of 
(p6[). The k th multicritical point is reached when t 2 k+i = — 2 k+i an< ^ a ^ °ther times are 
zero. In this case the string equation becomes 

VR k [u] = a k vx . (29) 

This is a 2k th order differential equation which as x — > oo has asymptotic solutions of the 
form 

oo 

v ~ x^(l + ^2vix~ l{2+ ^ ] ) , (30) 



which upon the identification k 2 = x ( 2+ s ) gives the genus expansion of the specific heat 

CO 

-(i + E/^ 23 )' (31) 



v 2 ~ X k 



h=l 



where f h = 2v h + J2i 1+ i 2=h v h v h- 

The connection to the r-function formalism of the mKdV hierarchy is shown by noting 
that the specific heat v 2 can be written in the form [T7]] 



v 2 = -d z log (nr 2 ) (32) 

with n and r 2 the r-functions of the mKdV hierarchy (pq) . These are simply connected to 
the Miura transformed functions wi = t> 2 +?/ and w 2 = f 2 — t/ by U{ = —2d 2 logr^ , i = 1, 2. 
Then the partition function is given by 

Z = n • r 2 , (33) 



which is the continuum analog of (|13|) . 

The Virasoro constraints [|17| are obtained by first substituting (p8|) into (|6[) and then 
using ( |32|) . The result is 

Lon = fin , (34) 

with L = ^2kL (k + |)t 2 fc+i at 2 fc + i an< ^ ^ an ar bitrary constant. The flows and the 

recursion relations relate L n+ i to L n and one obtains 

Ln7"i = with n > 1 , (35) 

where L n = V?° n (& + I)i2fe+ii57 — 1- h TVl—i 757 m ■ We will further discuss 

n <L^k=0\ 2> tn-ri- dt 2 ( k + n )+i 2 £—<K=\ dt 2 k-iot 2 ( n -k)+i 

the Virasoro constraints in section 5. 

4. The mKdV Hierarchy and the Sato Grassmannian 

As we already mentioned in the introduction, the analysis of the solutions of the 
string equation in the Sato Grassmannian Gr depends crucially on the association of the 
mKdV r-functions n and r 2 to points V\ and V2 in the big cell of the Sato Grassmannian 
Gr(°K In this section we take a pedestrian approach to explaining this association and the 
reader familiar with the subject might want to skip to the next section. For more rigorous 
treatments on the subject see |J2|| and the references therein. 



9 



Since the Sato Grassmannian is an infinite dimensional generalization of finite dimen- 
sional Grassmannians, we start by reviewing the relevant concepts in the finite dimensional 
case. For a nice review along these lines see | ]29|| . The Grassmannian Gr(k,N) consists 
of all /c-dimensional linear subspaces of C^. A point V G Gr(k,N) is described by a 
basis {vi} with i = 1, . . . , k and a basis of the orthogonal complement of V {uii} with 
i = k + 1, . . . , N. Then the pair (v,w) specifies a point in Gr(k,N). A pair (v',w'), 
however, gives the same point if 



(v r , w') = (v, w) 

Then 



A B 
C 



Gr(k,N) ~ GL(N)/P 

with P = Un Z) I The relation between Gr(k, N) and fennions is established by 



c t , 

considering the GL(iV) representation on a fermionic Fock space F defined by the vacua 

\k >= e± A . . . A e k < k\ = i ek . . .i ei < i\k >= 5 ik , (36) 

where {e^} is a basis of and i ei (ej) = Sij is the inner product operator. The fermionic 
operators are defined by 

ipl=e t A\ X > V] = * ei |x>, (37) 
and satisfy canonical anticommutation relations 

tyi, V]} = , {fa, ^} = {ijjj, = . (38) 
The vacua \k > carry charge k and tpli^i) create a charge +1(— 1). Then 

i/jl\k>=0 i = l,...,k ifj i \k>=0 i = k+l,...,N 
< k\ip\ = $ i = k + l,...,N <k\ipi = i = l,...,k, 

The Plucker embedding is defined by assigning to every point V € Gr(k, N) a state 

\v >= c Vx A . . . A v k with Vi = ■> (40) 

where {vi} is a basis of V and c is an arbitrary constant. A change of basis Vi — > ciijVj 
corresponds to c — > (deta) c and the state \v > is well defined. The condition 

^[ Vi ]\v>=0 Vi, (41) 



(39) 
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with ^[ui] = J2 v ij' l J J i defines equivalently the state \v > up to the constant c. 
Then a G gl(iV) acts on F by 

a\x>=J2^i ai ^ x> \X>^F, (42) 
and on the space of operators on F by 

[ipi, a] =^ a ik^k , [a, ipj] = ^2 ^l a ki ■ (43) 
k k 

The action of g G GL(iV) is defined by exponentiation of (|42| ). For example 

£ViV>I 2 • • -Mi 2 ■ • • |o >= (V^ktV^k ■ ■ ■ teVOiiteV'k • • • |o > (44) 

with (tp^g)i = tyjgji and (gVOi = Qij^j- Then a gl(iV) operator a acting on V £ Gr(k,N) 
by av = J2( a ij v j) e i corresponds to a fermionic operator a = ^'ipjaij^j- Then if a\ <-> ai 
and a 2 <-» 02, equations (^3|) give 

[ai,a 2 ] <-» [ai,a 2 ] • (45) 

Moreover note that if 

a|u >= const. \v ><^> aV C V. (46) 

The state |t> > belongs to the GL(iV) orbit of the state \k >. Since for \v >= v±A. . .Avk 
every vector V{ can be written in the form V{ = gei for some fixed g G GL(iV), we have 
that \v >= g\k > as defined in fl44T). Therefore the image of Gr(k,N) under the Plucker 
embedding can be identified with the orbit GL(iV)|/c >. 

The r-functions are given by fermion correlators 

r<? =< O > v =< k\0\v > , (47) 

with O a zero charge operator. Since the topology of Gr(k, N) is non-trivial, we divide 
it into cells (U a ,a G I). A point V G U a is represented by a basis {v\ } and the state 
\ v >(«)= v [ a) A...Auf. Then if V G W a n U b we have = a { k f ] vf ] and 

r^ a) =deta^rf 6) 

Therefore the r-functions are really sections of a determinant line bundle over Gr(k,N) 
whose transition functions are given by deta^ ab \ 
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Most of the results carry over almost unchanged to the infinite dimensional case. For 
the infinite dimensional vector space we consider the space of formal Laurent series 



H = {^2a n z n , a n = for n > } 

n 

and its decomposition 

H = H+ © H- , 

where H + = {Yl a nZ n , a n = for n 3> 0}. Then the big cell of the Sato Grass- 

n>0 

mannian Gr^ consists of all subspaces V C H comparable to H + , in the sense that the 
natural projection n + : V — * H + is an isomorphism. Then V admits a basis of the form 
{<pi(z)}i>o where <fii(z) = z l + lower order terms. The Plucker embedding ([|Dp is defined 
by the semi-infinite wedge product 

\v>=C(f) 1 (z) Afciz) A... . (48) 

Care has to be taken so that a GL(oo) change of basis <fii(z) — > aij(f)j(z) does not introduce 
infinities, since det a can be infinite. We choose a set of admissible bases for V G Gr to 
be those whose matrix relating {7r + (<^(z))}i>o to {z*}i>o differs from the identity by an 
operator of trace class. Then the fermionic representation is defined on the Fock space 
built on the vacuum state of zero charge 

|0>= 1 A z A z 2 A... , (49) 

by fermions tpj and tpi defined as in fl37|). The states (m > 0) 



|m>= V4--.'0i'|O>, I -ra>='0_ TO+ i...Vo|O> (50) 

are the filled states with charge m and — m respectively. The generalization of gl(iV) is 
given by gl(oo) and is represented on F by its central extension gl* (oo) with 

a = Yl : ^i a ^j : ( 51 ) 

where 

: 4^ ■■ = 4^3- < 4^ >= { *'f j , t 1 > < ° n (52) 



Ml i < 
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is the normal ordering. The reason for introducing normal ordering is that the naive 
operator J2i j i J i a ij' l J J j m aps an admissible basis to a non-admissible one. 

The connection of the fermion representation of Gr^ and the KP and mKP hierarchies 
is made explicit by making use of the boson- fermion equivalence in two dimensions. The 
fermionic currents 

Jn = ^ : ^n-r^r ■ U eZ (53) 

satisfy the bosonic commutation relations 

[Jm, Jn] = md m - n . (54) 

By representing the bosonic Fock space by B = C[t±,t2, ■ ■ ■ , ; u, the space of polyno- 
mials in ti,t2, . . • , ; u, it -1 , g|- and — nt_ n (with n > 0) act as creation and annihilation 
operators on B satisfying the algebra (|54l). Then fermionic operators can be mapped to 
operators acting on B and states in F to states in B by mapping the state \m > of F 
to u m . Then the k th modified KP hierarchy r-functions correspond to correlators (|4"7| ) 
where O = eYp>^ tpJp and the states \v > correspond to the GL(oo) orbit of \i > with 
i = 0, . . . , k — 1. In particular the solutions to the second mKP hierarchy is given by two 
r-functions 

n(t) =< i-l\expC^t p J p }g\i-l > (i = l,2), (55) 

where g G GL(oo). The modified KdV hierarchy that arises in UMM is the second reduced 
mKP hierarchy of the above equation and it corresponds to eliminating from the 
dependence on the even times {t2 n }- Therefore every solution T\(t) and r 2 (t) of the mKdV 
hierarchy corresponds to points V\{b) and V^(t) in GA ^ given by the states \vi(t) > = 
ex P{ Yl tpJp}9\i _ 1 >• Then the time dependence of Vi(t) is given by 

d 

Vi(t) >= J2k+i\vi(t) > and J2k\vi{t) >= 0, (56) 



dt2k+i 

or by using the correspondence 

d 



V l (t)=z 2k+1 V l (t) and z? k Vi{t)cVi{t). (57) 



dt2k+i 



Then Vi(t) = exp{J2t2k+iz 2k+1 }V l = 7 (t, z)V t . 
k 
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5. The Solutions to the String Equation 



Since to every solution of the mKdV hierarchy correspond points V±(t) and V 2 (t) in 
Q r {o) satisfying fl57|) , one would like to determine those that are solutions to the string 
equation (|25|). This is particularly easy because the commutator [P, Q-] is equal to a 
constant [p2| . 

Consider the space \1/ of pseudodifferential operators W = ^ Wi(x)d l where the 

i<k 

functions Wi(x) are taken to be formal power series (i.e. Wi(x) = Yl w ikX h , Wik = , k 

fc>0 

0). W is then a pseudodifferential operator of order k. It is called monic if Wk(x) = 1 
and normalized if Wk-i(x) = 0. The space \1/ forms an algebra. The space of monic, 
zeroth-order pseudodifferential operators forms a group Q. 
There is a natural action of \1/ on H defined by 

x m d n . jj ^ H 

Then it is well known |3(J that every point V G Gr^ can be uniquely represented in the 
form V = SH + with S G Q. This will imply that for every operator Q_ we can uniquely 
associate a pair of points V 1 ,V 2 e Gr(°). 

Indeed, consider Si and S2 G £? such that 

SQ.S" 1 = Q_ (58) 

where 

*-(? £)■«-- (2 0)' (59) 

Then 

S 1 {d + v)S^ 1 = d, 

(60) 

S 2 {d-v)Si 1 = d. 

Si and S2 can be shown to exist and are unique up to a redefinition 5 4 — > Si R with 
i? = 1 + Xli>o r i^~* an d n constants. 

since is given uniquely by V = SH +: the operator Q_ determines two 

spaces V\ = SiH + and V2 = S 2 H + . Conversely given spaces V\ and V 2 determine Q_ 

14 



uniquely. The operator <2_, however, is a differential operator and Vi, V 2 cannot be arbi- 
trary. Indeed, since every differential operator leaves H + invariant, we obtain 

(d + v)H+c H+ <*Si 1 dS 2 H+ c H+ 

V 2 c Vi (61) 

<^z V 2 C V\ 

Similarly, z V\ C V 2 . Notice that these conditions are consistent with the second equation 
in (|57D. 

The transformation SQ-S -1 = <2_ is a similarity transformation and the string equa- 
tion will be left invariant if we define Vn.) = SVr^S -1 . Then the solution to [P(k)> Q-] = 1 
is easily found to be given by 

fi(k) = f A q J , where A k = — + ^aiZ 21 and on =const. (62) 

The requirement that V be a differential operator is equivalent to the conditions 
Ak V\ C V 2 and Ak V% C V\. The space of solutions to the string equation is the space of 
operators Q- such that there exists V<k) with \P(k), Q-] = 1- We conclude that this space 
is isomorphic to the set of elements V\, V 2 C Gr^ that satisfy the conditions: 



(63) 



(64) 



zVtCV 2 zV 2 dVi 
A k V 1 cV 2 A k V 2 cV 1 

k 

for some Ak = 4- + ^ OLiZ 2% . 

i=0 

The string equation is left invariant by the flows ([57]). Indeed 

z70M)Vi C <y(t,z)V 2 ^zV 1 (t) c V 2 (t) 
A k {t) 7 (z, t)Vi C 7 (£, z)V 2 A k {t) Vi(t) C V 2 (t) , 

where 

A k {t) = -fAk-f' 1 = A k - J2( 2k + l)t2k+iz 2k (65) 

k 

and analogous equations with V\ and V 2 interchanged. 

It is now easy to see that implies the Virasoro constraints for the r-functions. 
Without going into the details (see [^3|), we first notice that the operators l n = z 2n+1 A 
leave Vi invariant 

z 2n+l AViCVi. (66) 
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Then using the correspondence (|5), one can construct the corresponding fermion opera- 
tors l n and from them their bosonic counterparts L n |[23|| . These have the exact form as 
equations (^4|)-(|35|). We can immediately see that they form a Virasoro algebra by noting 



that L ~ z 2n+1 f- are the generators of the Virasoro algebra and by using lemma ( f45[) . 



Since l n leave Vi invariant, then using ([|6[) we conclude that the operators L n annihilate 
the r-functions t\ and t<i and obtain equations (|33J) and . 

We conclude this section by showing how conditions ( |5"5D can be used to calculate the 
space of solutions to the string equation [P2[ . We will start by describing the spaces V±,V2- 



First choose vectors fa(z) , fa(z) G V\, such that 

fa(z) = 1 + lower order terms , fa(z) = z + lower order terms . 
Then the condition z 2 V\ C V\ and tt + (Vi) = H + shows that we can choose a basis for V\ 

(f)l, fa, Z 2 fa, Z 2 fa, . - - 

Since z V\ C V2 and 7r + (V2) = H + we can choose a basis for V2 to be 

ijj, Z(f)i, Z(f) 2 , Z 3 (j)i,Z 3 (j)2, ■ ■ ■ 

where ip(z) = 1 + lower order terms. Using z V2 C V\ we have zij) = a<pi + (3 cf>2- Choose 
(pi, fa such that zip = fa. Then we obtain the following basis for Vi, V2 {(j> = fa): 

Vi : 0, zip, z 2 (p, z s ip, . . . 

(67) 

V 2 : tp, zfa z 2 tp, z 3 (p, . . . 

Then it is clear that (p,ip specify the spaces Vi,V2. Using the conditions AV\ C V2 and 
AV2 C V\ we obtain 

(4- + fk(z 2 ))4> = Poo(z)(p + Poi(z)iP 

dz (68) 

(— + f k (z 2 ))iP = P 10 (z)<P + P n (z)iP . 
dz 



Since a generic system of the form ( p8[) will lead to exponential evolution of the 
functions (p and ip, the requirement that they keep their polynomial form puts severe 
conditions on Pij(z). A detailed calculation shows that the space of solutions to the string 

equation (^5]) is the two fold covering of the space of matrices ^Pij(z)^ with polynomial 

entries in z such that Pq\{z) and Pio(z) are even polynomials having equal degree and 
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leading terms and Poo (z) and Pu(z) are odd polynomials satisfying the conditions Pqq(z) + 
Pu(z) = and degP Q0 (z) < degP 01 (z). 
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